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We measure the topological susceptibility by cooling with an over-improved action. In contrast with usual 
cooling, large instantons survive over-improved cooling indefinitely. By varying the parameter of the over-improved 
cooling action, we measure the instanton size distribution. 



1. Motivation 

A variety of methods are available to measure 
the topological charge Q of a non-abelian lattice 
gauge field. To choose among them, one is guided 
by generally conflicting goals: 

1) the measurement of Q should be insensitive to 
"dislocations" , short-range fluctuations caused by 
the lattice discretization; 

2) Q represents a topological index, and thus 
should be an integer; 

3) the measurement of Q must be "objective" : it 
must not depend on a tuning parameter whose 
value is fixed by subjective criteria. 

The most common methods, "geometric" , "an- 
alytic" and "cooling", miss at least goals 1, 2, 
and 3 respectively. Here we try to address the 
deficiencies of the latter. 

Cooling (or smearing) consists of an iterative, 
local minimization of the Wilson action. Cool- 
ing dynamics can be slowed down if desired, by 
under-relaxation [1]. Cooling has become widely 
used to gain qualitative insight into the large- 
scale features of the gauge field. Under cool- 
ing dynamics, many observables show metasta- 
bility: their value reaches a plateau before the 
configuration eventually collapses to the trivial 
vacuum. How soon this plateau is reached and 
how long it lasts depends, however, on the observ- 
able, and even more problematically, on the initial 
gauge configuration. This lack of robustness has 
prevented widespread acceptance of cooling as a 
quantitative method. 



The behavior of instantons under cooling is well- 
known. As the lower curve of Fig.l shows, the 
action of a lattice instanton decreases monoton- 
ically with its size. Thus in the process of mini- 
mizing the action, cooling will shrink instantons 
until their size becomes 0(a): then the lattice 
discretization is so coarse that the non-trivial 
topology of the gauge-field cannot be detected, 
even by a geometric method. The instanton has 
"dropped through the lattice". When one looks 
for a plateau in the topological charge Q under 
cooling, one tries to observe the slow shrinking of 
the largest instantons. It is not at all clear how 
many smaller, physical instantons have already 
disappeared at this stage. 
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Figure 1. Action of a lattice instanton versus its 
size, for the Wilson action (bottom curve) and 
various over-improved actions S(e). 
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To remedy this situation we propose to use over- 
improved (01) cooling [2]. An over- improved ac- 
tion, at tree-level, makes positive the first term 
0((a/p) 2 ) in the instanton action: so the action 
of large instantons now increases as they shrink. 
Since by construction the lattice instanton action 
must go to zero with its size, the action will de- 
velop a maximum for a size p pea k- The upper 
curves of Fig.l show this effect for various values 
of the over-improvement parameter e. Minimiz- 
ing this action by cooling will then shrink instan- 
tons of size smaller than p pea k until they disap- 
pear; those larger than p pea k will expand to the 
maximum size allowed by the finite lattice. They 
will be stable against any number of cooling steps. 
Thus we can consider an objective procedure to 
measure the net topological charge contributed by 
instantons of size p > p pea k'- cool with an over- 
improved action until a stable state is reached; 
measure Q then. This measurement is facilitated 
by the fact that instantons are now maximally 
large and smooth: the simplest, analytic method 
(measuring J F^F^) can be employed. By vary- 
ing the over-improvement parameter e, and thus 
Ppeak(e), the size distribution of instantons can 
be investigated. 

2. Results 

For economic reasons, we test our approach on 
SU(2) rather than SU(S). We use the action [2]: 

S(e) = ^Tr(l-W lxl )+ € -^Tr(l-W 2x2 )(l) 

where e < for over- improvement. This action 
has been used so far to produce instantons in situ- 
ations where they exist in the continuum, ie. with 
twisted boundary conditions. Here we want to 
use it with periodic b.c, where no continuum in- 
stanton exists. We were therefore concerned that 
our stable state after 01 cooling might be dis- 
torted from a continuum instanton. Fig. 2 shows 
the distribution of the topological charge den- 
sity away from the instanton center, measured by 
Fp.vF^v, where F M „ is extracted from the 01 ac- 
tion (1), on a 16 4 lattice after (infinitely) long 
cooling with e = — 1. The solid line is t'Hooft's 
continuum ansatz oc ( x 2 P +p 2 ) 4 [periodic images 
have a negligible effect on the ansatz]. The fit is 
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Figure 2. Profile of the topological charge density 
for the maximum size instanton (L = 16, e = —1). 
The solid line is t'Hooft's ansatz with p/a = 3.5. 

excellent up to large distances where the bound- 
ary of the box is felt. The total charge and total 
action are very close to their continuum values 
(the maximum deviation was ~ 7% on the small- 
est 8 4 lattice). The surprise however is that the 
instanton size does not exceed p/a ~ 3.5, even on 
a 16 4 lattice. We consistently observed a max- 
imal size Pmax/a < L/4 in all our simulations. 
This can be contrasted with much larger sizes 
observed with twisted b.c. (see eg. [3]). Thus 
periodic b.c. induce very strong finite size effects 
on instantons. 

We kept periodic b.c. nonetheless, to compare 
our results with previous work [1]. Our procedure 
consisted then of generating a sample of config- 
urations with the Wilson action, and of cooling 
each of them several times, with 01 actions of 
successively more negative e. If n(p) is the num- 
ber of instantons of size p, we measure on each 
configuration, after (infinitely) long 01 cooling: 

Q(e) = / dp n(p) (2) 

The upper bound ~ L/4 in the integral comes 
from the finite volume as just discussed. The 
lower bound is the threshold size for instanton 
survival under 01 cooling: it is p pea k (e) as visible 
on Fig.l. Then from the distribution of Q(e) over 
our sample, we extracted an effective susceptibil- 
ity 

XL(e) = < Q 2 (e) > / V (3) 
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Figure 3. Topological susceptibility as a function 
of e, for various lattice sizes ({3 = 2.4). 

Our method is expensive: we wanted to make 
sure that we reached a stable state, so we cooled 
with 100-700 sweeps for each value of e. This is 
the price we pay for an objective method. 

Our sample consisted of 1000 (540) configu- 
rations for L = 8, 10, (12) at = 2.4, sepa- 
rated by 50 Monte Carlo sweeps. Remarkably, 
no measurable autocorrelation could be detected 
in (2(e), even for e = —0.1: the largest instan- 
tons appeared as drawn from a random sequence. 
Our MC update was a 9-to-l mixture of over- 
relaxation and heatbath. Preliminary tests with 
a pure heatbath update revealed no autocorre- 
lation either. This is a welcome finding on the 
efficiency of pure gauge MC updates, in sharp 
contrast with the long autocorrelation times ob- 
served for Q with dynamical fermions, using Hy- 
brid Monte Carlo [4]. 

Our results for the effective susceptibility (3) are 
shown in Fig. 3. Several observations are in order: 
1) as e is changed, the lower cutoff p m m(e) is 
changed. The effect on the susceptibility is dra- 
matic (a factor ~ 10 on the 8 4 lattice). If we 
could push Pmin to zero, presumably we would 
recover the susceptibility obtained via the geo- 
metric method. It is clear that measurements of 
the topological susceptibility are meaningless if 
not accompanied by an estimate of the lower cut- 



off Pmin • 

Ideally p m in should be adjusted to include all 
physical instantons but leave out all lattice arti- 
facts. Although an estimate p/a ~ 1.93 for this 
ideal value can be obtained [3], it really should 
be extracted a posteriori, by looking at the scal- 
ing with (3 of the instanton size distribution: lat- 
tice artifacts will give a divergent contribution be- 
low Pmin , while the distribution of real instantons 
above p m i n will scale. 

2) as L is changed, the upper cutoff p m ax(L) is 
changed, more or less in proportion to L. There- 
fore the difference \L\ ( e ) — Xl 2 (e) should be in- 
dependent of e: it comes from large instantons, 
which fit inside the large box but not the small 
one. Fig. 3 shows clearly that this is the case: the 
data for L = 8 lie below the rest, with a nearly 
constant offset. On the other hand there is lit- 
tle difference between L = 10 and 12, indicating 
that the contribution of large instantons of size 
> Pmoa;(10) ~ 2.3 to the true, infinite volume 
susceptibility is very small. 

3) using a _1 (/3 = 2.4) = IJGeV, an arrow 
marks the phenomenologically plausible suscep- 
tibility (200MeV) 4 . Such a result could be ob- 
tained with e ~ —0.3, which is the value favored 
in [5]. 

The next step consists of mapping our results 
as a function of p m i n , not e. p m in{t) is the av- 
erage minimum size of an instanton, generated 
by the Wilson action, necessary for survival un- 
der 01(e) cooling. Here we measured p m i n {t) by 
monitoring 5(e) evaluated on a shrinking instan- 
ton, being cooled with the Wilson action. The 
size was determined in a manner similar to [1], by 
looking for a hypersphere containing half the to- 
tal charge. One difficulty is that an instanton for 
the Wilson action is not an instanton for the OI 
action (it does not satisfy the classical equations 
of motion) : so different Wilson instantons of the 
same size may suffer different fates under 01(e) 
cooling. In addition, p m j„ also depends on the 
proximity of the instanton to a lattice site. An 
attempt at averaging is performed in [3]. Here 
we merely acknowledge these effects as a source 
of uncertainty. 

Extracting p TO m(e) from Fig.l, the effective sus- 
ceptibility is displayed in Fig. 4 as a function 
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causes a bias, which depletes the small size pop- 
ulation and may be responsible for the observed 
peak. Our procedure avoids this bias. 



Figure 4. Effective susceptibility versus minimum 
instanton size. 

of Pmin- The value (200MeV) 4 corresponds to 
Pmin/a ~ 1.96, very close to the estimated cutoff 
for lattice artifacts ~ 1.93 [3]. We also present 
cubic fits (higher degree polynomials are not con- 
strained enough): the maximal instanton sizes 
given by these fits (~ 2.25 and ~ 2.35) are con- 
sistent with the sizes actually measured after 01 
cooling. 

Finally we can extract the instanton size distri- 
bution. Since Q 2 (e) = J dpidp2n(pi)n(fi2), and 
assuming that instantons interact only weakly, 
so that correlations between n(pi) and n{pi) are 
negligible, we get 

<<9 2 (e) >= I Pmin (e) d P< n2 (p) >. 
so that 



< n 2 (p min ) > 



dpr> 



<Q 2 >. 



The differentiation of our cubic fit (L = 12) is 
shown in Fig. 5. It is striking that, although al- 
lowed by the degree of the polynomial fit, and 
unlike all previous studies, no maximum appears 
in the instanton size distribution. Several expla- 
nations are possible: 

• we have large statistical uncertainties, espe- 
cially in the determination of p m i n (e) • 

• our assumption < n(p\)n(p2) >= S{pi, Pi) may 
be flawed because of (anti-)instanton interactions, 
non-negligible in such a small box. 

• all other studies measure instanton sizes on 
cooled configurations. At that stage, small in- 
stantons have already shrunk and disappeared, 
and small, close-by A-I pairs have annihilated. 
Measuring the size distribution at that stage 
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Figure 5. Fitted instanton size distribution 
< n 2 (p) > (L = 12). 
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